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1. Introduction 



The top quark is the heaviest fundamental particle discovered so far. Its interactions are still 
relatively unexplored and their experimental investigation may lead to exciting new discover- 
ies. To mention only a few examples: due to its large mass, it is not excluded that the top quark 
decays into yet unobserved particles like charged Higgs bosons or supersymmetric particles - 
or, if these particles are heavier than the top quark, they may mediate top quark decay, leading 
also to new decay modes and/or branching ratios that differ from the Standard Model (SM) pre- 
dictions (for overviews, see, e.g., refs. EI3)- Specifically, the V— A structure of the top quark 
decay vertex is modified in many extensions of the Standard Model BHD. In the production of 
top quarks, the resonant production of heavy spin-O-particles could lead to interesting signatures 
J5J|SJ|7|. Its large mass makes the top quark also a good probe of the electroweak symmetry 
breaking mechanism. In this context it will be important to check whether the Yukawa coupling 
of the top quark is as predicted by the Standard Model. Finally, the question whether the discrete 
symmetry CP is violated in top quark production and decay has been the subject of numerous 
investigations (cf. ref. for a review). The foreseen large data samples at the upgraded Teva- 
tron (~ 10 3 — 10 4 ff/year) and at the LHC (~ 10 7 /f/year) will tell us more about these issues in 
the near future. Of help in this context will be another special feature of the top quark, namely 
that the spin properties of top quark pairs are predictable by perturbation theory. This is due 
to its large decay width, Tf M ~ 1.5 GeV ^> Aqcd, which serves as a cut-off for hadronization 
effects. In particular, the top quark decays so fast that the information about its polarization is 
not diluted by hadronization but transferred to the decay products. Thus observables related to 
the spins of the t and t quarks can be constructed and used reliably for the detailed study of 
the dynamics of top quark production and decay [9|. Apart from searching for non-standard 
effects, the study of the t and t spin properties is interesting even within the framework of the 
SM: They probe the 'quasi-free' nature of the top quark, thus allowing us to study properties 
of a 'bare' quark. Further, a measurement of spin correlations would provide a lower bound on 
\V t t,\ without assuming the existence of three quark generations [10]. 

Here we confine ourselves to top quark pair production and decay at hadron colliders and inves- 
tigate the top quark polarization and spin correlation phenomena that are induced by the strong 
interaction dynamics. We discuss predictions for the normal polarisation of the top quarks and 
tt spin correlations at next-to-leading order (NLO) in the QCD coupling a s . 



2. Theoretical framework 

Needless to say, for a correct interpretation of upcoming and future data on top quark produc- 
tion and decay at the Tevatron and the LHC, precise theoretical predictions within the SM are 
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needed. We consider here the following reactions: 
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h x h 2 ^tt + X^{ B _ J ' , (1) 



where = p,p; £ = e,jji,x, and j t (jf) denotes jets originating from hadronic t (t) decays. 
Experimental analysis of the above processes requires predictions of the fully differential cross 
sections. The calculation of these cross sections at next-to-leading order QCD simplifies enor- 
mously in the leading pole approximation (LPA), which amounts to expanding the full ampli- 
tudes for the reactions listed in eq. ([T]) around the complex poles of the top quark propaga- 
tors. Only the leading pole terms are kept in this expansion, i.e., one neglects terms of order 
T t /m t ~ 1%. Within the LPA, the radiative corrections can be classified into factorizable and 
non-factorizable contributions. Here we consider only the factorizable corrections; for the non- 
factorizable contributions see ref. ifTTl . Further we apply the on-shell aproximation for the top 
quark propagators: 

r/m^o k — m A + iml ml 

The necessary ingredients at NLO QCD within this approximation are the differential cross 
sections for the following parton processes 

qq — > tt, gg — > tt (to order oc^), (3) 

qq -> tig, gg -> tig (to order a^), (4) 

gq(q) — > tiq(q) (to order 0$), (5) 

t^b£v,bqq' (to order a s ), (6) 

where we have to keep the full information on the t and i spins. 



3. Spin density matrix 

The fully differential cross section J0f ac t. for the production and subsequent decay of top quark 
pairs at the parton level can be written in terms of a spin density matrix R and decay density 
matrices p, p: 

dGfact. = Tr»,f spins {R P p) • (V) 

The unnormalized spin density matrix is explicitly given by the following expression (i = 

qq,gg, •••): 

tfaM'P' = E('(^ v)i{k h p)x|<r \i) (t(k t , a')i(k h &)x\*r\i)*, (8) 
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where the sum runs over all unobserved degrees of freedom. The decomposition of the spin 
density matrix with respect to the t and t spin spaces reads: 

R=Al®l + B + o<g>l + l®oB +Qj a ! '®a J '. (9) 

The polarization of the top quark (antiquark) is encoded in B ± , e.g. 

B+_TrNl] 
P f _2(S,)- — - Tt[r] , (10) 

where S t denotes the top quark spin operator. The matrix C encodes the spin correlations of the 
top quark and antiquark: 

Ms,A,j)= l Tr[R] ' J =f ■ (ii) 



3.1. Normal polarization 

If only strong interactions are taken into account then the polarization of t and t in pp.pp — > ttX 
can only be normal to the event plane due to parity invariance of QCD. Normal polarization 
requires absorptive parts in the scattering amplitude, i.e. one-loop diagrams with discontinuities. 
For the two initial states i = qq.gg, 

Pj = P' F = ^(y, 1 ?)A J (12) 

where n is the unit vector normal to the event plane and y = cos with 9 denoting the scattering 
angle of the top quark in the parton center-of-mass frame. The normal polarization at the parton 
level is a percent effect lfT2l fT3l . Several observables to study this effect were proposed and 
studied in ref. [fT2l. At the Tevatron, and probably even at the LHC, the normal polarization of 
the top quark induced by QCD will be very difficult to observe. This means that it provides a 
sensitive probe of non-standard strong rescattering effects in hadronic top quark pair production. 



3.2. Spin correlations at leading order 

The correlation between two observables 0\, O2 is defined as 

with dOi = yj (Of) — (Oj) 2 . Assuming parity invariance we have for hadronic production of top 
quark pairs at leading order (i.e., no absorptive parts): 

(s;> = (5/) = (i4) 
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and therefore 



Q 



con(S' n Si)=4(S' t Sl} = ^-. 



For the process qq — > tt, the spin correlation matrix at LO is quite simple: 
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yyk t + y/l-y 2 K 



(15) 



(16) 



(17) 



where [3 = (1 — Am 1 / s) l l 2 , J = (1 — |3 2 )~ 1//2 , and k t denotes the direction of the top quark 
Further, k = p g — yk and d = p g — (p ? • d)d, where p q denotes the quark direction. 

For qq — > ?F the direction d is the optimal spin basis at leading order QCD, because at LO: 



diCffdj 
Aw 



1. 



(18) 



For any [3 and y, the tt spins are 100% correlated w.r.t. to this basis ff4l. Following the nomen- 
clature of ref. ifbUl we shall call this basis the off-diagonal basis in the following. For the 
Tevatron, where this choice of spin axis is useful there is an equally efficient but simpler possi- 
bility. At threshold, the top quark pair is in a 3 5i state and we have 



d >Vq- 



(19) 



This suggests that at the Tevatron the direction p of, say, the proton beam is an equally good 
choice of spin axis. 



At high energies, helicity is conserved and we have 

-3^1- 
d P -^k t . 



(20) 



For the process gg — > tt, no optimal spin basis exists. The LO expression for the spin correlation 
matrix is quite lengthy, and we therefore discuss only the limiting cases. At threshold, the top 
quark pair is in a 1 Sq state and 
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--Sy. 



At high energies, we have 



rss „ , 

Cy 3^1 1 „ 2 



Ass 3 lJ 1 \ +y 2 



and helicity conservation is reflected by 
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K t,i K t,j ^ 'J 



(21) 



(22) 



(23) 
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4. Observing spin correlations 



The tt spin correlations show up in certain angular distributions/correlations of the top decay 
products, e.g. for the dileptonic decay channel t — > £ + vb, t — ► £'~vb the following distribution 
is sensitive to the correlations: 

1 d 2 a(hih 2 -> ttX -> e+i'-X) 1 , 

^-V — ^ ^ = T (l+^icose + +B 2 cose_-Ccose + cose_). (24) 

a acos0 + cos0_ 4 

In eq. (1241 . 9+,9 are the angles of £ ± in the ? (F) rest frame with respect to arbitrary spin 
quantization axes a, b, e.g.: 

a = — b = k f (helicity basis) , 

a = b = p (beam basis), (25) 
a = b = d (off-diagonal basis) . 

The coefficient C reflects the strength of the tt spin correlations for the chosen quantization axes, 

OCD OCD a ^ 

— 1 < C < + 1 1 9 1 . Further we have = Bf" = if a, b are chosen to be in the production 
plane due to P invariance of QCD. 



5. Spin correlations at NLO 

Within the LPA, the coefficient C of eq. d2"4"l) factorizes: 

C = k + k_D (26) 

with the tt double spin asymmetry 

n = a(TT)+a(ll)-g(TI)-a(IT) m , 

o(n)+o(u)+o(u)+oar)' 

In eq. (|2"71) . a(|T) denotes the hadron cross section for top quark pairs with t(t) spin parallel to 
the chosen spin quantization axis a(b) etc. The numbers K± are the spin analysing powers of 
charged leptons in decays t(t) — > b(b)£ ± v{v). The decay distribution reads 

^y^V = ^ (1 + ±K±cos^ ± ) , (28) 
rJcosf>± 2 

where -&± are the angles of £ ± w.r.t. the t (F) spin. 

5.1. Spin analysing power of top quark decay products 

If the t or t quark decays hadronically, one can use other decay products as spin analysers. One 
defines in analogy to eq. (l28t for t — > bW + — > b£ + v or bqq': 

ra = ^ 1+K / C0S *)' <29) 
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The leading order results for K/ are given in Table 1 . In order to compute the spin correlation 



Table 1 : Leading order results for the spin analysing power of top quark decay products. In the 
last column, 'qq' jet' stands for the least energetic non-b-quark jet in hadronic t decays hi 51/ . 





£+, d, s 


V, u 
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W+ 




K / 


1 


— 0.31 


-0.41 


0.41 


0.51 



coefficient C in NLO in a?, we need the QCD corrections to the spin analysing power K/. For 
leptonic decays, the corrections are tiny I fTol 

K + = K_ = l-0.015a„ (30) 

implying that the charged lepton is a perfect analyser of top quark spin. However, since only 
about 5% of the decays of tt are purely leptonic (e,//), it is also important to compute the QCD 
corrections for hadronic decays t — > bqq' . The results are H7] (using (X s (m t ) = 0.108): 

Kj = 1- 0.57a, = 0.94, (31) 
K b = -0.41 x (1 -0.34a,) = -0.39, (32) 
Kj = +0.51 x (1 -0.65a,) = +0.47, (33) 

where K,- is the analysing power of least energetic non-fr-quark jet. One observes that the d, s- 
jets have the highest analysing power, but the reconstruction of their direction is very difficult 
due to the low efficiency of charm tagging. A better choice is to use the b-jet or the least 
energetic non-^-quark jet. The NLO results for C in the single lepton channel are less sensitive 
due to the factor Kj j,, but this loss in sensitivity is overcompensated by higher statistics, since 
one has about 30% (e + /j) +X single lepton /f decays. (Nonleptonic top quark decays at order 
a, were also analyzed in ref. [18].) 



5.2. Double spin asymmetries at NLO at the parton level 

The spin dependent cross sections that enter the double spin asymmetry (l27t are calculated as a 
convolution which reads schematically: 

o(TT) =PDF / s®6(TT),... (34) 

We renormalize the top mass in the on-shell scheme, and a, in the modified minimal subtraction 
(MS) scheme. Factorization is performed in the MS scheme, and we set /jp = /jr= /j. The 
results at NLO QCD for the MS subtracted parton cross sections qq — > tt(g), gg — > tt(g), and 
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q(q)g — > ttq(q) with tt spins summed over was computed already more than 10 years ago lfT9l 
|20l|2T]|. It can be written as follows: 



a(s,m 



&(TT)+o(U) + a(U)+6(|t) 

/ (0) (T1) +4mJ/ 1 )(ri) +/ 1 )(Ti)ln( A / 2 /m f 2 ) 



a 2 



nr, 



(35) 



where T| = §/ (4m 2 ) — 1. In our calculation of the spin-dependent cross sections refs. Il22ll23l 
we obtain these cross sections as special cases and we find perfect agreement with the results of 
refs. lfT9ll20ll2T1l . An analogous decomposition can be defined for the numerator of the double 
spin asymmetry: 



6D = 6(TT) + 6(U)-6(U)-6(!T) 

{g(°)(Tl) +47KX, [gW(Tl) +g( 1 )(T])ln( i u 2 /m 2 )] } 



a] 



(36) 



The following figures display our results for the functions g^°\g^ and g^ calculated for sev- 
eral spin quantization axes. In eq. (l36l) the coupling a,, denotes the six-flavour coupling of - . 
Fig. 1 shows our results for the helicity basis. In Figs. 2 and 3 corresponding results for the 
beam basis and the off-diagonal basis are displayed. 
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Figure 1: Left: Scaling functions g(°> (t|) (dotted), g^^(T|) (full), and g^ l \r\) (dashed) in the 
helicity basis for the process qq — ► tt(g). Middle: The same for the process gg — > tt(g). Right: 
The functions g^\r\) (full), and g^\y\) (dashed) for the process qg — > qtt It22\ 12317 . 



It is natural to express the above partonic cross sections in terms of the MS coupling a{ in 
/ = 6 flavour QCD. However, for the evaluation of hadronic observables, e.g. eq. (l34b . the 
parameter change 0Cj — > a s ~ using the standard MS relation 
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Figure 2: The same as Fig. 1, but for the beam basis l[22\\23\l . 
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Figure 3: The same as Fig. 1, but for the optimal basis KT2\ l23l/ . 



a{ 6 M =oc{ 5 (^) 



(37) 



is necessary, in order to make contact with the physics and formalism incorporated in the PDF 

lig 



libraries. The evolution of of 5 is determined by the beta function with nl ght = 5 flavours. 
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5.3. NLO results for differential decay distributions 



We now turn to the coefficient C in the double angular distribution d24l . For definiteness we 
discuss the spin correlation C only for the dilepton channels. These results were obtained in ref. 
[ 24 1 . The numbers in Table 2 are obtained for /jf = pr — m t =175 GeV and using the PDF's 
CTEQ5L (LO) and CTEQ5M (NLO) E3. At the Tevatron, the dilepton spin correlations are 

Table 2: Coefficient C of the double distribution to leading and next-to-leading order in 
CL s for the helicity basis, the beam basis (where the proton beam is taken as the spin quantiza- 
tion axis) and the off-diagonal basis. The parton distribution functions of ref. H25V were used 
choosing the renormalization scale jjr equal to the factorization scale /up = m t = 175 GeV. 





pp at y/s = 2 TeV 


pp at y/s = 14 TeV 


LO NLO 


LO NLO 


Chel. 
Cbeam 
C ff. 


-0.456 -0.389 
0.910 0.806 
0.918 0.813 


0.305 0.311 
-0.005 -0.072 
-0.027 -0.089 



large in the beam and the off-diagonal basis. Thus one sees there is practically no difference 
between these two choices as far as the sensitivity to QCD-induced spin correlations is con- 
cerned; yet the beam basis may be simpler to implement in the analysis of experimental data. 
The QCD corrections are about —10%. At the LHC, the beam and off-diagonal bases are bad 
choices (due to the dominance of gg — > tt). Here the helicity basis is a good choice, and the 
QCD corrections are small. The inclusion of the QCD corrections reduces the dependence of 
the tt cross section on the renormalization and factorization scales significantly. The same is 
true for the product aC, as can be seen from Fig. 4. Table 3 shows the dependence of the NLO 
results for C on the scale p (upper part) and on the choice of the PDFs (lower part). At the 
Tevatron the spread of results for different PDFs is larger than the scale uncertainty: The results 
for C using the CTEQ5 and MRST98 distributions agree up to a few percent, but the difference 
between GRV98 and MRST98 at the Tevatron is about 10%. The main reason for this strong 
dependence on the PDFs is that the contributions from the gg and the qq initial state enter with 
a different sign. This offers the interesting possibility to constrain the PDFs by measuring tt 
spin correlations. 

In all results above we used m t = 175 GeV. A variation of m t from 170 to 180 GeV changes 
the results for the Tevatron, again for p = m t and PDFs of ref. E51 as follows: Ch e i. varies from 
-0.378 to -0.397, C bea m from 0.790 to 0.817, and C off . from 0.797 to 0.822. For the LHC, 
Chei. changes by less than a percent. 
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Figure 4: Left: Dependence o/oCbeam at LO ( dashed line ) and at NLO ( solid line ) on [i = pr — 
Hf for pp collisions at y/s = 2 TeV, with PDF of ref. / l25l/ . Right: Same, but for aChei. for pp 
collisions at y/s = 14 TeV. 



Table 3: Upper part: Dependence on the scale /j of the correlation coefficients C, computed 
at NLO with the PDF of ref.Kl5\l. Lower part: Correlation coefficients C at NLO for p. = m t 
and different sets of parton distribution functions: GRV98126], CTEQ5[25], and MRST98 (c- 
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The results above have been obtained without imposing any kinematic cuts. Standard cuts on 
the top quark transverse momentum and rapidity only have a small effect on C: For the Tevatron, 
demanding |k^ | > 15 GeV and \r t j\ < 2 leads to the following results: Ch e \. = —0.386, Cb ea m = 
0.815, C ff. = 0.823. For the LHC, when imposing the cuts |k^ ? | > 20 GeV and \r t j\ < 3, we 
find Chei. =0.295. 
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A tool for simulating these spin correlations at leading order in the QCD coupling exists ll28l . 
Apart from the double differential distribution (|24b there are also other distributions which are 
sensitive to spin correlations [29|. Finally let us mention that there has been a first measurement 
of spin correlations in the off-diagonal basis by the DO collaboration OOl . It is based on six 
dilepton events from Run I. They find 

C ff. > -0.25 @ 68% confidence level (38) 

This demonstrates that top quark spin correlations can be studied already at the Tevatron. It is 
expected that the correlations can be established at the 2a level using Run II data. 



6. Conclusions 



QCD-induced spin correlations of top quark pairs produced at hadron colliders are large effects. 
They can be studied at the Tevatron and LHC. For the Tevatron the QCD corrections to the 
leading order predictions are sizeable but under control. The degree of correlation depends on 
the choice of the spin quantization axis. We have shown that using, at the Tevatron, the direction 
of one of the hadron beams as spin quantization axis is as efficient as the off-diagonal axis which 
has received much attention in the literature. Spin correlations are suited to study in detail 
the interactions of top quarks. As we have pointed out they should, in particular, be a useful 
tool for constraining PDFs. Taking the PDFs, once they have been determined with sufficient 
precision, as input, spin correlations will be an important tool for the search for new effects in 
top quark pair production. Future work on the theory side will include the implementation of 
the NLO matrix elements in an event generator, a study of non-factorizable corrections, and the 
resummation of Sudakov-type large logarithms for the spin-dependent cross sections. 



Acknowledgments 



A.B. would like to thank the organizers of the Epiphany conference for their hospitality. 



References 

[1] M. Beneke et al., (2OOO), |hep^ph/0003033l 

[2] ECFA/DESY LC Physics Working Group, E. Accomando et al., Phys. Rept. 299 (1998) 1, 

|hep-ph/97054421 

[3] A. Soni and R.M. Xu, Phys. Rev. Lett. 69 (1992) 33, 

[4] PL. Cho and M. Misiak, Phys. Rev. D49 (1994) 5894, |hep-ph/9310332| 



11 



[5] W. Bernreuther and A. Brandenburg, Phys. Rev. D49 (1994) 4481, |hep-ph/ 93122101 
[6] D. Dicus, A. Stange and S. Willenbrock, Phys. Lett. B333 (1994) 126, |hep-ph/9404359| 

[7] W. Bernreuther, M. Flesch and P. Haberl, Phys. Rev. D58 (1998) 1 14031, 

|hep-ph/97092841 

[8] D. Atwood et al., Phys. Rept. 347 (2001) l, |hep-ph/0006032[ 
[9] J.H. Kiihn, Nucl. Phys. B237 (1984) 77, 

[10] T. Stelzer and S. Willenbrock, Phys. Lett. B374 (1996) 169, |hep-ph/ 9512292[ 

[11] W. Beenakker, F.A. Berends and A.P Chapovsky, Phys. Lett. B454 (1999) 129, 

|hep-ph/99023041 

[12] W. Bernreuther, A. Brandenburg and P. Uwer, 
Phys. Lett. B368 (1996) 153, |hep-ph/9510300| 

[13] W.G.D. Dharmaratna and G.R. Goldstein, Phys. Rev. D41 (1990) 1731, 

[14] G. Mahlon and S. Parke, Phys. Lett. B411 (1997) 173, |hep-ph/9706304[ 

[15] A. Czarnecki and M. Jezabek, Nucl. Phys. B427 (1994) 3, |hep-ph/9402326[ 

[16] A. Czarnecki, M. Jezabek and J.H. Kiihn, Nucl. Phys. B351 (1991) 70, 

[17] A. Brandenburg, Z.G. Si and P. Uwer, Phys. Lett. B539 (2002) 235, |hep-ph/02050231 

[18] M. Fischer et al., Phys. Rev. D65 (2002) 054036, |hep-ph / 1 1 3 2~2l 

[19] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303 (1988) 607, 

[20] W. Beenakker et al., Phys. Rev. D40 (1989) 54, 

[21] W. Beenakker et al., Nucl. Phys. B351 (1991) 507, 

[22] W. Bernreuther, A. Brandenburg and Z.G. Si, Phys. Lett. B483 (2000) 99, 

|hep-ph/0004184[ 




[23] W. Bernreuther, A. Brandenburg, Z.G. Si and P. Uwer, Phys. Lett. B509 (2001) 53, 

|hep-ph/01040961 

[24] W. Bernreuther, A. Brandenburg, Z.G. Si and P. Uwer, Phys. Rev. Lett. 87 (2001) 242002, 

|hep-ph/01070861 

[25] CTEQ, H.L. Lai et al., Eur. Phys. J. C12 (2000) 375, |hep-ph/9903282[ 

[26] M. Gliick, E. Reya and A. Vogt, Eur. Phys. J. C5 (1998) 461, |hep-ph/9806404[ 




12 



[27] A.D. Martin et al., Eur. Phys. J. C4 (1998) 463, |hep-ph/ 9803445L 

[28] S.R. Slabospitsky and L. Sonnenschein, Comput. Phys. Commun. 148 (2002) 87, 

|hep-ph/0201292[ 

[29] W. Bernreuther, A.Brandenburg, Z.G. Si and P. Uwer, in preparation. 
[30] DO, B. Abbott et al., Phys. Rev. Lett. 85 (2000) 256, |hep-ex/0002058[ 



13 



